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We discuss exclusive double diffractive (EDD) production of Higgs boson and heavy quark - heavy 
antiquark pairs at high energies. Differential distributions for cc at \fs = 1.96 GeV and for bb at 
= 14 TeV are shown and discussed. Irreducible leading-order bb background to Higgs production 
is calculated in several kinematical variables. The signal-to-background ratio is shown and several 
improvements are suggested by imposing cuts on b (6) transverse momenta and rapidities. 
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I. INTRODUCTION 



pair production pp — > p(qq)p in the color singlet state as 



Exclusive production of the Higgs boson can be an al- 
ternative to the present studies of Higgs in inclusive pro- 
cesses. There is recently a growing theoretical interest in 
studying exclusive processes. Only a few processes have 
been measured so far, mostly at the Tevatron (see [l| and 
references therein). Khoze, Martin and Ryskin developed 
an approach in the language of off-diagonal unintegrated 
gluon distributions. This approach was applied to exclu- 
sive production of Higgs boson [2j . In our recent papers 
we applied the same formalism to exclusive production 
of cc and bb quarks. Quite large cross sections have been 
found @||. 

The cross section for the Standard Model Higgs pro- 
duction is of the order of 1 fb for light Higgs [2J . The dom- 
inant bb decay channel is therefore preferential from the 
point of view of statistics. It was argued that the leading- 
order contribution is rather small using a so-called J z — 
rule. Here we show a quantitative calculation which 
goes beyond this simple rule. In our calculation we in- 
clude exact matrix element for massive quarks and the 2 
— > 4 phase space. This fully four-body calculation allows 
to impose cuts on kinematical variables. Different types 
of backgrounds to Higgs production were studied before 
e.g. in Ref.p. 



II. FORMALISM 

Let us concentrate on the simplest case of the produc- 
tion of qq pair in the color singlet state. Color octet state 
would demand an emission of an extra gluon which con- 
siderably complicates the calculations. We do not con- 
sider the qqg contribution as it is higher order compared 
to the one considered here. 

We write the amplitude of the exclusive diffractive qq 
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where \ q , \q are helicities of heavy q and q, respectively. 
Above fi S and f 2 s are the off-diagonal unintegrated 
gluon distributions in nucleon 1 and 2, respectively. 

The longitudinal momentum fractions of active glu- 
ons are calculated based on kinematical variables of 
outgoing quark and antiquark: x\ = exp(+t/ 3 ) + 

^ exp(+7/ 4 ) and x 2 = ^ exp(-y 3 ) + ^exp(-y 4 ), 
where 7713,4 and m 4j t are transverse masses of the quark 
and antiquark, respectively, and 7/3 and 7/4 are corre- 
sponding rapidities. 

The bare amplitude above is subjected to absorption 
corrections. The absorption corrections are taken here in 
a simple multiplicative form. 

Let us consider the subprocess amplitude for the qq 
pair production via off-shell gluon-gluon fusion. The ver- 
tex factor V^ 2 _ — V^ 2 _(qi,q 2 , fci, k 2 ) in expression (|2.ip 
is the production amplitude of a pair of massive quark 
q and antiquark q with helicities X q , Xq and momenta 
fci, k 2 , respectively. The color singlet qq pair production 
amplitude can be written as Q 

V \ q l\ (9i . 92 , fci , k 2 ) = 77+77; V^l 2 9 ' ^ (qi ,q 2 ,k 1 ,k 2 ), 

The tensorial part of the amplitude reads: 
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The coupling constants g s — > g s {^1 i)<7s(A*r 2)- I n the 
present calculation we take the renormalization scale to 
be p? r l = p? r 2 = Mgg/4 or M^. The exact matrix ele- 
ment is calculated numerically. Analytical formulae are 
shown explicitly in Q. 
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The off-diagonal parton distributions (i=l,2) are cal- 
culated as 

fi ^( x ii Qi,t> A 1 >^i) = 

^9 dl^itf !*?=<% * W 

where 51/2(9?, A* 2 ) is a Sudakov-like form factor relevant 
for the case under consideration. It is reasonable to take 
a running (factorization) scale as: p\ = /i 2 . = or 

M 2 - 

The factor i? 9 here cannot be calculated from first prin- 
ciples in the most general case of off-diagonal UGDFs. 
It can be estimated in the case of off-diagonal collinear 
PDFs when x' <C x and xg = x~ A (l — x) n . Then 

R 9 = 2 ^ 3r r?A+4f ■ Typically R g ~ 1.3 - 1.4 at 
the Tevatron energy. The off-diagonal form factors are 
parametrized here as F(t) — exp(_B fft). In practical 
calculations we take B h = 2 GeV~ 2 . In the original 
KMR approach the following prescription for the effective 
transverse momentum is taken: Q\ t = min (g 2 t , g 2 t ) 
and Q\ t = min (g 2 t , q\ t ). In evaluating f-y and f% 
needed for calculating the amplitude (|2.ip we use differ- 
ent collinear distributions. It was proposed to express 
the S1/2 form factors in Eq. (|2.3I) through the standard 
Sudakov form factors as: 

S 1/2 (qln 2 ) = ^T g (qln 2 )- (2-4) 

The cross section for the four-body reaction is calcu- 
lated as 

da = 5j|A^2^4| 2 (27r) 4 (5 4 (p a + p b - Pl - p 2 - p 3 - Pa) 

ri 3 pi d 3 p 2 d 3 p 3 d 3 p 4 

(2tt) 3 2,Ei (27r) :i 2B 2 (2tt) 3 2_E 3 (27r) b 2,E4 - 

The details how to conveniently reduce the number of 
kinematical integration variables are given elsewhere. 

III. RESULTS 

A. pp — > ppec 

Let us proceed now with the presentation of differen- 
tial distributions of charm quarks produced in the EDD 
mechanism. In this case we have fixed the scale of the 
Sudakov form factor to be p, — M cS /2. Such a choice 
of the scale leads to a strong damping of the cases with 
large rapidity gaps between q and q. 

In the left panel of Fig. [I] we show distribution in ra- 
pidity. The results obtained with the KMR method are 
shown together with inclusive gluon-gluon contribution. 
The effect of absorption leads to a damping of the cross 
section by an energy-dependent factor. For the Tevatron 
this factor is about 0.1. If the extra factor is taken into 
account the EDD contribution is of the order of 1% of 
the dominant gluon-gluon fusion contribution. 

The corresponding rapidity-integrated cross section at 
y/s = 1960 GeV is: 6.6 pb for the exact formula, 2.4 pb 



for the simplified formula (see Eq. (12.3[1 ). For comparison 
the inclusive cross section (gluon-gluon component only) 
is 807 ph. 
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FIG. 1: Rapidity distribution of c or c (left) and transverse 
momentum distribution of c or c (right). The top curve is for 
inclusive production || while the two lower lines are for the 
EDD mechanism with leading-order collinear gluon distribu- 
tion [101 ]. The solid line is calculated from the exact formula 
and the dashed line for the simplified formula 0]. An extra 
cut on the momenta in the loop Qt,cut = 0.26 GeV 2 was im- 
posed. Absorption effects were included by multiplying the 
cross section by the gap survival factor Sg = 0.1. 

In the right panel of Fig. Q] we show the differen- 
tial cross section in transverse momentum of the charm 
quark. Compared to the inclusive case, the exclusive con- 
tribution falls significantly faster with transverse momen- 
tum than in the inclusive case. 



B. pp — > ppbb 

In parallel to the exclusive bb production, we calculate 
the differential cross sections for exclusive Higgs boson 
production. Compared to the standard KMR approach 
here we calculate the amplitude with the hard subprocess 
9*9* ~~ ^ H taking into account off-shellness of the active 
gluons. The details of the off-shell matrix element can 
be found in Ref. Q- In contrast to the exclusive Xc pro- 
duction [8J, due to a large factorization scale ~ Mh the 
off-shell effects for g*g* — > H give only a few percents. 

The same unintegrated gluon distributions based on 
the collinear distributions are used for the Higgs and 
continuum bb production. In the case of exclusive Higgs 
production we calculate the four-dimensional distribu- 
tion in the standard kinematical variables: y, i 1 , £2 and 
<t>. Assuming the full coverage for outgoing protons we 
construct the two-dimensional distributions da/dycPpt in 
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FIG. 2: The bb invariant mass distribution for y/s = 14 TeV 
and for -2.5 < y < 2.5 corresponding to the ATLAS/CMS 
detectors. The absorption effects were taken into account by 
multiplying by the gap survival factor So = 0.03. 
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FIG. 3: The bb invariant mass distribution for t/s = 14 TeV 
and for b and b jets in the rapidity interval —2.5 < y < 2.5 
corresponding to the ATLAS detector. The absorption effects 
for the Higgs boson and the background were taken into ac- 
count by multiplying by the gap survival factor Sg = 0.03. 
The left panel shows purely theoretical predictions, while the 
right panel includes experimental effects due to experimental 
uncertainty in invariant mass measurement. 



Higgs rapidity and transverse momentum. The distribu- 
tion is used then in a simple Monte Carlo code which 
includes the Higgs boson decay into the bb channel. It 
is checked subsequently whether b and b enter into the 
pseudorapidity region spanned by the central detector. 

In Fig.[5]we show the most essential distribution in the 
invariant mass of the centrally produced bb pair, which is 
also being the missing mass of the two outgoing protons. 



In this calculation we have taken into account typical de- 
tector limitations in rapidity —2.5 < yb,yi < 2.5. We 
show results with different collinear gluon distributions 
from the literature: GRV [TO], CTEQ [U, GJR [II] and 
MSTW [ll]. The results obtained with radiatively gen- 
erated gluon distributions (GRV, GJR) allow to use low 
values of Q t — got > Qit , lit whereas for other gluon dis- 
tributions an upper cut on Q t is necessary. The lowest 
curve in Fig J2] represents the 77 contribution [4j. While 
the integrated over phase space 77 contribution is rather 
small, it is significant compared to the double-diffractive 
component at large M bb > 100 GcV. This can be under- 
stood by a damping of the double diffractive component 
at large M bb by the Sudakov form factor 0, H[ . In ad- 
dition, in contrast to the double-diffractive component 
the absorption for the 77 component is very small and 
in practice can be neglected. 

In the left panel of Fig(3]we show the double diffractive 
contribution for a selected (CTEQ6 11]) collinear gluon 
distribution and the contribution from the decay of the 
Higgs boson including natural decay width calculated as 
in Ref. 
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see the sharp peak at M bb = 120 GeV. The 
phase space integrated cross section for the Higgs pro- 
duction, including absorption effects with Sq — 0.03 is 
somewhat less than 1 fb. The result shown in Fig J3] in- 
cludes also the branching fraction for BK(H — > bb) s» 0.8 
and the rapidity restrictions. The second much broader 
Breit-Wigner type peak corresponds to the exclusive pro- 
duction of the Z° boson with the cross section calculated 
as in Ref. [15j. The exclusive cross section for y/s = 14 
TeV is 16.61 fb including absorption. The branching frac- 
tion BR(Z° — > bb) s» 0.15 has been included in addition. 
In contrast to the Higgs case the absorption effects for 
the Z° production are much smaller [ll| . The sharp peak 
corresponding to the Higgs boson clearly sticks above the 
background. In the above calculations we have assumed 
an ideal measurement. 

In reality the situation is, however, much worse as both 
protons and in particular b and b jets are measured with a 
certain precision which automatically leads to a smearing 
in M b i . Experimentally instead of M b i one will measure 
rather two-proton missing mass (M pp ). The experimen- 
tal effects are included in the simplest way by a convo- 
lution of the theoretical distributions with the Gaussian 
smearing function G(M) = -y= — exp ^ ^ M ~^*l H ^ j with 

(j = 2 GeV which is determined mainly by the preci- 
sion of measuring forward protons. In the right panel we 
show the two-proton missing mass distribution when the 
smearing is included. Now the bump corresponding to 
the Higgs boson is below the bb background. With the 
experimental resolution assumed above the identification 
of the Standard Model Higgs seems rather difficult. The 
situation for some scenarios beyond the Standard Model 
may be better. 

Can the situation be improved by imposing further 
cuts? In Fig. 2] we show the distribution for the EDD 
background in y b and y b . In contrast to the Higgs [J] 
the cross section for the bb continuum has maxima far 



4 



EDD pp— >ppbb 1 ! 7.6 <^ B < 122.4 GeV 




from the diagonal. This can be used to impose cuts on 
quark/antiquark rapidities. In Fig.[5](left panel) we show 
the result for a more limited range of b and b rapidity, 
i.e. not making use of the whole coverage of the main 
LHC detectors. Here we omit the Z° contribution and 
concentrate solely on the Higgs signal. Now the signal- 
to-background ratio is somewhat improved. This would 
be obviously at the expense of a deteriorated statistics. 
Similar improvements of the signal-to-background ratio 
can be obtained by imposing cuts on jet transverse mo- 
menta. Detailed studies of the role of cuts is discussed in 



FIG. 4: Two-dimensional distribution in yi, and yi for the 
EDD bb background. 
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FIG. 5: The bb invariant mass distribution for \/s = 14 TeV 
for a limited range of b and b rapidities: — 1 < y < 1. 



I am indebted to Rafal Maciula and Roman Pasechnik 
for collaboration on the issues presented here. 



[1] M.G. Albrow, T.D. Coughlin and J.R. Forshaw, 

Arxiv.1006.1289. 
[2] V. A. Khoze, A. D. Martin and M. G. Ryskin, Phys. Lett. 

B 401, 330 (1997); 

A. B. Kaidalov, V. A. Khoze, A. D. Martin and 

M. G. Ryskin, Eur. Phys. J. C 33, 261 (2004). 
[3] R. Maciula, R. Pasechnik and A. Szczurek, Phys. Lett. 

B 685, 165 (2010). 
[4] R. Maciula, R. Pasechnik and 

A. Szczurek, arXiv: 1006.3007 [hep-ph], in print in 

Phys. Rev. D. 

[5] R. Maciula, R. Pasechnik and 

A. Szczurek larXiv: 101 1 .58421 [hep-ph] . 
[6] M. Chaichian, P. Hoyer, K. Huitu, V.A. Khoze and AD. 

Pilkington, JHEP 0905 (2009) 011. 
[7] R. S. Pasechnik, O. V. Teryaev and A. Szczurek, Eur. 



Phys. J. C 47, 429 (2006). 
[8] R. S. Pasechnik, A. Szczurek and O. V. Teryaev, Phys. 
Rev. D 78, 014007 (2008); 

R. S. Pasechnik, A. Szczurek and O. V. Teryaev, 
[9] T.D. Coughlin and J.R. Forshaw, JHEP 1001, 121 (2010). 
[10] M. Gliick, E. Reya and A. Vogt, Z. Phys. C67, 433 
(1995). 

[11] J. Pumplin et al., JHEP 0207, 012 (2002). 

[12] M. Gliick, D. Jimenez-Delgado, E. Reya, Eur. Phys. J. 

C53, 355 (2008). 
[13] AD. Martin et al., Eur. Phys. J. C63, 189 (2009). 
[14] G. Passarino, Nucl. Phys. B 488, 3 (1997). 
[15] A. Cisek, W. Schafer and A. Szczurek, Phys. Rev. D80 

074013 (2009). 



